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Abstract 

The problem of controlling surge and stall in jet engine compressors is of fundamental importance 
in preventing damage and lengthening the life of these components. In this paper, we use the Moore- 
Greitzer mathematical model to develop an output feedback controller for these two instabilities (only 
one of the three states is measurable). This problem is particularly challenging since the system is 
not completely observable and, hence, none of the output feedback control techniques found in the 
literature can be applied to recover the performance of a full state feedback controller. However, we 
show how to successfully solve it by using a novel output feedback approach for the stabilization of 
general stabilizable and incompletely observable systems. 

1 Introduction and Problem Description 

In this paper we consider the problem of controlling two instabilities which occur in jet engine compressors, 
namely rotating stall and surge. Rotating stall develops when there is a region of stagnant flow rotating 
around the circumference of the compressor causing undesired vibrations in the blades and reduced 
pressure rise of the compressor. Surge is an axisymmetric oscillation of the flow through the compressor 
that can cause undesired vibrations in other components of the compression system and damage to the 
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engine. In |l8| , Moore and Greitzer developed a three-state finite dimensional Galerkin approximation 
of a nonlinear PDE model describing the compression system. Since its development, several researchers 
have used the Moore-Greitzer three state model (MG3) to design stabilizing controllers for stall and 
surge. The available control approaches may be divided into three main categories: 1) Linearization and 
linear perturbation models (e.g., ^ ^ among others); 2) Bifurcation analysis (e.g., ^, [l^, |l|); 

and 3) Lyapunov based methods (e.g., |3[ |2^). Most existing results focus on the development of state 
feedback controllers, thus complicating their practical implementation as in [Q, where the authors use 
sensor arrays (2D sensing) to implement a state feedback control law depending on the squared amplitude 
of the first harmonic of asymmetric flow and the derivative of the air flow through the compressor. In 
a partial state feedback controller simplifies practical implementation by only requiring measurements 
of the mass flow and plenum pressure rise (hence 2D sensing is not needed). On the other hand, the 
limitation of this partial state feedback controller lies in the fact that it cannot globally stabilize a unique 
equilibrium point. 

To the best of our knowledge, no attempt has been made to design a stabilizing output feedback 
controller (using only plenum pressure rise feedback) based on a full-state feedback control law. This is 
probably due to the fact that MG3 becomes unobservable when there is no mass flow through the com- 
pressor, i.e., the system is not uniformly completely observable (UCO), and none of the techniques found 
in the output feedback control literature (e.g., |2l|, |2^, |^, ^ |l6|, ||) can be employed for the solution 
of this problem. In this paper we introduce a new globally stabilizing full state feedback control law for 
MG3, and we employ the theory developed in jlj, |l^ for the output feedback control of incompletely 
observable nonlinear systems to regulate stall and surge by using only pressure measurements. The MG3 
model is described by (see ^ for an analogous exposition) 

$ = + «'c($) - 3$i? 



vi/ = ^($-a>T) (1) 



R = aR{l - $2 _ ^(0) > 
where $ represents the mass flow, is the plenum pressure rise, i? > is the normalized stall cell 
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squared amplitude, $t is the mass flow through the throttle, cr = 7, and /? = \/\/2. The functions ^'c($) 
and $t(^) are the compressor and throttle characteristics, respectively, and are defined as ^^ci^) = 
*Co + 1 + 3/2$ - f/2$^, * = i(f + where is a constant and 7 is the throttle opening, 

the control input. Given the static relationship existing between $7^ and 7, without loss of generality, 
in what follows we will design a controller assuming that <^t is our control input. Our control objective 
is to stabilize system (|l|) around the critical equilibrium = 0,$'^ = 1,5*^ = ^fci^"^) — 'I' Co + 2, 
which achieves the peak operation on the compressor characteristic. We shift the origin to the desired 
equilibrium with the change of variables (/) — ^ — l,tp ~ '9 — ^^Co ^ 2. System then becomes 

-aR^ -aR{2(j) + (l3^) 
^ = -il;-3/2(f>^ ~l/2(f>^ -3R(j)-3R (2) 

The pressure rise (and hence ■0) is the only measurable state variable. It is readily seen that this system 
is input output feedback linearizable with relative degree one (the first derivative of ■0 contains the input 
7), and its zero-dynamics are nonminimum phase. 



2 State Feedback Control Design 

For convenience, in the remainder of the paper we will redefine the control input to be u = $t — 1- 
Next, notice that Assumption A2 in jl^, |l2j is satisfied since, for example, a stabilizing control law for 
(^ is given in |^ by means of backstepping design. However, the control law proposed in ^ turns 
out to be quite complex. In it is shown that a linear partial state feedback control law of the type 
u — diip — d2(f> achieves either a unique asymptotically stable equilibrium point with domain of attraction 
{(i?, 0, ^) e M"^|i? > 0} or two equilibria on the axisymmetric and stall characteristic, with domains of 
attraction {(i?, t^,?/;) e R^\R = 0} and {(i?, t^,?/;) e R^|i? > 0}, respectively (see Theorem 3.1 in g). 
Here, this problem is overcome by viewing system (|^) as an interconnection of two subsystems, namely the 
i?-subsystem and the {(p, ■!/')-subsystem, and then building a full state feedback controller which makes the 
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origin of (||) an asymptotically stable equilibrium point with domain of attraction {(i?, (/>, i/;) G M'^ji? > 0}, 
as seen in the next theorem. 

Theorem 1 For system with the choice of the control law 

{1- (3'^kik2)^ + P^k2ijj + 3(3^kiR(j) (3) 
where ki and k2 are positive scalars satisfying the inequalities, 
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k,>^+'—^ (4) 

L-(C(T + 3)2>0 (5) 
9,2 9fci (kl - 1)2 



f/ie origin is an asymptotically stable equilibrium point with domain of attraction A = {{R, <p, ip) S M^|i? > 
0}. 

Proof. For the sake of simplicity, redefine the control input to be u' = — ^ (u — 0) , so that the last 
equation in becomes ijj = u' . Next, notice that system (H) can be viewed as the interconnection of 
two subsystems: 



[Si] R = -cTi?^ [S: 



A Lyapunov function for [Si], defined on the domain {R e M | i? > 0}, is Vi = i?, and its time derivative is 
readily found to be Vi = —aR^ thus showing that the origin of [Si] is an asymptotically stable equilibrium 
point of [5*1], and its domain of attraction is {R e R | i? > 0}. As for subsystem [S2] the analysis found in 
Section 2.4.3 in Q suggests using V2 = + ^4''^ + ^(0— kiipY , where fci is a positive design constant. 
Furthermore, in B, a stabilizing control law for [^2] is found to be u' = —cicj) + C2'4', where ci and 



4 



C2 are two appropriate positive constants. In the following we will show that, in order to stabilize the 
interconnection of systems [Si] and [82], one needs to add to u' — —ci4>+C2ip a term which is proportional 
to the product R(j). Based on these considerations, consider the following candidate Lyapunov function 
for system (||), 

V^CVi + V2^ CR+]-<f + .^04 _^ i (V; - fci0)' (8) 
Z a A 

where C > is a scalar. After noticing that V is positive definite on the domain A, and letting 
ip = ip — ki(j), we calculate the time derivative of V as follows, 



V = - CaR^ - CaR{2(j} + 0^) + ( + ^(f^) ( -V' - f 0^ - \'f - 3R(j) - 3R 

(9) 

+ + kiij) + ^kicf^ + ^^kicjr' + ikiRcj) + 3fcii? 

Here, as in we use the identity — — \4>^ = —5 + 4> + %4> to eliminate the potentially 
destabilizing term —(</) + A:i/20'^) 3/2(/>^. Next, substituting (||) into (^) (after taking in account the 
definition of u'), letting ki — ki — 9/8, and using the definition of tp, we get 



+ V' ( -(^2 - ki)i! + kl<f) + ^ki(j)^ + + UiR 



(10) 



Now notice that the expression — (0 + ■^f/''^) 5 (0 + |) can be discarded since it is negative definite, and 
that the term ^(jy^'ij) cancels out. After collecting the remaining terms, we get 

V<- CaR^ - {2Ca + i)R4> - [Ca + i)R4>'' - - + ^r] 0^ - ^i?^^ + 

^ ^ (11) 

+ i> (-{k2 - fci)V^ + (k\ - + ^fci02 + 3fcii? 
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By using Young's inequality five times we have 



3fci^^3 ^3hfR^^ , 



4 4 z 4 4A;i 



(/c^ - 1)0^ < (/.^ + 



Applying the inequalities above to (ni]) we get 



V < 



< 



R 



Ca 



2 2 
3 



4 1 4A:i 4 1 ^ 



i(Ca + 3-|fcil 



i (Ca + 3- ffci 



R 



^^"^^-4^^-iir 



(12) 



Hence, is negative definite on the domain A, provided that the quadratic form above is positive definite 
and that the coefficients multiplying 0^ and tp'^ be positive. By imposing the positive definiteness of the 
quadratic form we obtain Ca — | > 0, (Ccr — |) jkiki — j (Ca + 3 — > 0, while by imposing 

the positivity of the coefficients of the remaining two terms we get fci > (^'^'^+'^) — |- k2 > ki + jkf + 
||i _|_ ihi^D__ By using the definition of ki, inequalities (Q), (||), (||), and (|^) follow. In conclusion, if ki, 
^2, and C are chosen so that (^)-(0) hold, we have that V is negative definite on A which contains the 
origin. This leads to the conclusion that {R = 0, (/> = 0, '0 = 0} is an asymptotically stable equilibrium 
point, which in turn implies that {R = 0, = 0, t/i = 0} is an asymptotically stable equilibrium point. 
Our next objective is to show that ,4 is a region of attraction for the origin. This, however, is not 
immediately evident from our result, since the set e M.^\V < K,K > 0} is unbounded 

and, due to the presence of the term CR in V, it is not completely contained in A. In other words, it 
may happen that, while the Lyapunov function is decreasing, R becomes negative, and thus the state 
trajectory exits the set A, where V is guaranteed to be negative definite. Therefore, in order to complete 
our analysis, we need to show that A is invariant, which, together with < 0, implies that the set 
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{[R, (pjtpy G M.^ \ V < K, K > 0} n A is a region of attraction of the origin for any K > 0. This is readily 
seen by noticing that, on the boundary of i? = 0. From (||), R = imphes R — 0, thus proving that 
no trajectory of the system can cross the boundary of A, and therefore A is invariant. In conclusion, 
given any initial condition [i?(0), (/)(0), ?/'(0)]^ in A, there exists a constant A' > such that the initial 
condition is contained in the set {[R,ip,tp]^ e M'^ | y < K,K > 0} A, thus proving that the origin of 
system is an asymptotically stable equilibrium point with domain of attraction A. 




Figure 1: Comparison between the partial state feedback controller developed in [|| and the full state 
feedback controller (0). 



Remark 1: By using inequalities (§)-(0), it is easy to show that the only equilibrium point of the 
closed- loop system on the set A is the origin, as predicted by Theorem 1. Figure |l] shows the evolution of 
the closed-loop trajectories under the partial state feedback controller developed in |^ and the controller 
(||) for a particular choice of the coefficients di , ^2 , fci , A:2 . The partial state feedback controller stabilizes 
an equilibrium point different from the origin {R,4>,ip) = (0,0,0). 



Remark 2: Inequalities (^-(@) represent conservative bounds on ki and k2- In practical implementa- 
tion, these parameters may be chosen significantly smaller after some tuning. 
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In order to complete the state feedback design, we have to add an appropriate number of integrators 
at the input side of the system (see FoUowing the procedure outhned in we form the 

observabihty mapping Ti. 



(-U - ^ - 3 /202 _ l/2</)3 - 3i?0 - 3i?) 



(13) 



Notice that the observabihty assumption Al in ||TJ, is satisfied, for aU (j) 7^ —1, with = 2 in that 
given ye,M, and u, one can uniquely find R,(j),^p. The operating point = —1 corresponds to $ = 0, 
i.e., no mass flow through the compressor which is a condition we would like to avoid during normal 
engine operation. Since n„ = 2, we extend the system with two integrators ii — Z2, Z2 — v, u — zi. To 
simplify the notation in the following, define x — [R, (f>, ip]'^ , and rewrite (^) as i = f{x) + g{x)zi. Next, 
we find a stabilizing control law for the extended system by using the integrator backstepping lemma: 
V = Q!-zi-/c4Z2 = ip{x,z), where zi = zi-u, a = -kszi- ^g{x) + ^[f{x)+g{x) zi], Z2 = Z2-a, and 
k^, fc4 are arbitrary positive constants. This completes the design of a stabilizing state feedback for the 
extended system. The Lyapunov function of the closed- loop extended system isV = + Notice 

that, following the same reasoning as in the proof of Theorem 1 the set {[R, cj), ip, zi, Z2\^ e | i? > 0} 
is invariant; hence by applying the backstepping lemma we guarantee that the origin of the extended 
system is asymptotically stable with domain of attraction P = ^ x M^. 



3 Output Feedback Design 

The validity of the observability assumption Al in ||lj, |l^ allows us to design a stable observer. As already 
pointed out, Assumption Al in [|l^ |l^ is satisfied on the domain Xy.U — {[i?, (/>, -0] G K'^ | (/> > — l} x M^. 
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We first design the observer developed in , 



R = -aR^ - aR{2(j)- 



ih/p) + + 3i? + (3/2)<^2)(Z2/p2) + f3\l^/p^) 



3(1 + ^) 

= -V^ - 3/2 02 - 1/2 03 _ - 3i? + P^{l2/p'^)i^ - i^) 

Zl- (f) 



(14) 



/32 



{h/p)(^-^P) 



where p is a positive design parameter and the vector L ~ [hjh, h]^ G R'^ is chosen to be Hurwitz. Next, 
we calculate the solution P of the Lyapunov equation P{Ac — LCc) + [Ac — LCc)^ P = — /, where (Ac, Cc) 
is a canonical observable pair. In order to confine the observer estimates to within the observable space, 
we implement the following projection. 



m 

dx 



dn 



', N{o(^N{i,zyi + N,{izyz) 

V{ilz,z)^l^ NiizyVNiiz) 



if N{i z^i + N,{i, zYz > and I e dC^{z) 



otherwise 



where r = {S£')-^{S£')-^ , S ^ S'^ denotes the matrix square root of P, | = 7i:(.t, z), | = {^^+ ^^}' 
and C^{z) is the cube 



Q(z)=|eeM3|eie [ai,6i],6e 



-^(zi + 02), --^(21-62) 



^(-22-a3),-^(-Z2 + 63) 



which, when 02 < 1, is contained in TL{X , z), for all z (the scalars a^, hi, i = 1, 2, 3 have to be chosen to 



satisfy Assumption A3 in |14, IJ]). Finally, N{£^, z) and iVz(Cj z) are the normal vectors to the boundary 
of C^{z) with respect to ^ and z, respectively, and are given by 



[1,0, 0]T if 6 =61 [-1,0, 0]T if 6 =ai 

[0, 1, 0]T if I2 = -^(^1 - b2) [0, -1, 0]T if 6 = 'jjizi + 02) 

[0,0, 1]T if 6 = ■^i-Z2 + b3) [0,0, -1]T if ^3 = ^(-22-03) 
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0,0]T if =61 or Ci =ai 

^,0]^ ife2--/^(^i-fo2) 



^,0 



if 6 = ~7fT(2i + a2) 
if 6 = -0t{-Z2 - 03) 



Thus, the output feedback controller design is completed by letting v = (p{x^,z), and Theorem 2 in 
|l^ guarantees that the origin of the closed- loop system , controlled by w, is asymptotically stable 



A 



with domain of attraction V x fl^^ , where fl^^ ~ { [i?, -0] \V < C2, and i? > 0}, C2 > is the largest 
scalar such that Sl^^ C {[R, (p, ip]^ G M"^ | > —1}, and V C fl^^ can be made arbitrarily close to H.^^ by 



choosing p in (llj) small enough (see Theorem 2 in (Til E2j). 



4 Simulation Results 

Here we present the simulation results when the output feedback controller developed in the previous 
section is applied to system (^. We choose fci — 25 and ^2 = 1.1 • 10'"^ to fulfill inequalities in 
Theorem 1. In order to choose the size of the compact set C^{z) so that Assumption A3 in |l4| |l^ is 
satisfied, we may use the Lyapunov function V to calculate ft^^, choose C2 small enough to guarantee 
that C X, and use H to calculate bounds on ^ when x e fi^a- However, a more practical way to 
address the design of C^{z) consists of running a number of simulations for the closed- loop system under 
state feedback corresponding to several initial conditions [R.{Q),(j>{Q),tp{Q)]^ , and calculating upper and 
lower bounds for tp, (p, and —-0 — 3/20^ — l/2(f>'^ — SRcf) — 3R: these will provide the values of ai,bi, 
i — 1,2,3, respectively. By doing that, we found that whenever [i?(0), 0(0), -0(0)]^ G fio = {[R,4':i^V ^ 
E^li? G [O,O.1],0 e [-0.1, 0.1], -0 e [-0.5,0.5]}, we have that = -2, bi = 1, aa = -0.5, 62 = 1, 
03 — —0.5, 63 — 0.3 satisfy Assumption A3 in |l4[ |l^. We must point out that our choice of Uq is rather 
conservative and is made primarily for the sake of illustration. The actual domain of attraction V under 
output feedback control is larger that flo. In Figure |2| system and controller states, together with the 
control input, are plotted for two decreasing values of p confirming the theoretical predictions about the 
arbitrary fast rate of convergence of the observer found in Theorem 1 in |14]. Furthermore, the figures 
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R and its estimate 



(|> and its estimate 



V|/ and its estimate 




p = 0.1. 



R and its estimate 



(|i and its estimate \|/ and its estimate 

0.05 



Control input, Integrator state 

■ ■ 1 400 I ■ 




/9 = 0.02. 



Figure 2: Output feedback trajectories. 



also show the operation of the projection which prevents the observer from peaking and guarantees that 
(j) > —0.5. FinaUy, note that the output feedback trajectories approach the state feedback ones, as showed 
in Figure |[ 
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